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ABSTRACT 

The results from weak gravitational lensing analyses are subject to a cosmic variance error 
term that has previously been estimated assuming Gaussian statistics. In this letter we address 
the issue of estimating cosmic variance errors for weak lensing surveys in the non-Gaussian 
regime. 

Using standard cold dark matter model ray-tracing simulations characterized by f2 m 
()..'). S1 A = 0.7, h = 0.7, erg = 1.0 for different survey redshifts z s , we determine the 
variance of the two-point shear correlation function measured across 64 independent lines of 
sight. We compare the measured variance to the variance expected from a random Gaussian 
field and derive a redshift-dependent non-Gaussian calibration relation. 
We find that the ratio between the non-Gaussian and Gaussian variance at 1 arcminute can be 
as high as ~ 30 for a survey with source redshift z s ~ 0.5 and ~ 10 for z s ~ 1. The transition 
scale d c above which the ratio is consistent with unity, is found to be i) c ~ 20 arcmin for 
z s ~ 0.5 and d c ~ 10 arcmin for z s ~ 1. We provide fitting formula to our results permitting 
the estimation of non-Gaussian cosmic variance errors and discuss the impact on current and 
future surveys. 

A more extensive set of simulations will however be required to investigate the dependence 
of our results on cosmology, specifically on the amplitude of clustering . 
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1 INTRODUCTION 

Weak lensing by large scale structure, i.e. cosmic shear, offers a 
direct way of investigating the statistical properties of matter in 
the Universe, without making any assumptions on the relation be- 
tween dark and luminous matter. Current surveys are large enough 
to provide high precision constraints on cosmology and the lat- 
est measurements performed with the Canada France Hawaii Tele- 
scope Legacy Survey (Hoekstra et al. 2006, Semboloni et al. 2006) 
is a step in that direction. Most of the cosmological constraints 
from weak lensing use two-point shear statistics i Refregie r 2003 1 
I Van W aerbeke & Mellier 2003 1, and a crucial step in these cos- 
mological parameter measurements is the estimate of error bars 
and systematics. Several papers address, statistically, the is- 
sue of systematics from E and B modes (Crittenden et al. 200U 
IPen et al. 20021 |Schneider & Kilbinger 2006) , but only few pa- 
pers address the estimation of cosmic variance of cos- 
mic shear measurements (White & Hu 2000) |Cooray & Hu 200T| 
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Schn eider et al. 20 02 ). The latter assumes that the error on the two- 
point shear correlation function follows Gaussian statistics. How- 
ever, we know that this is not the case at small scales where non- 
linear effects become important. Cooray & Hu 2001 use the dark 
matter halo model in Fourier space to study non-Gaussian covari- 
ance. A tentative calibration of this effect on the aperture mass 
statistic ( Van Waerbeke et al. 2002 ) showed that departure from 
Gaussianity is expected to occur at angular scales 10 arcmin- 
utes. The purpose of this Letter is to estimate the non-linear co- 
variance of the two-point shear correlation function in real space, 
such that it can be of direct practical use for weak lensing studies, 
as in Schneider et al. 2002, without having to calculate high or- 
der correlation functions semi-analytically. Using ray-tracing sim- 
ulations for a model close to the concordance cosmological model 
( Spergel et al. 20061 at different source redshift slices, we obtain a 
redshift dependent calibration formula of the Gaussian covariance 
derived in Schneider et al. 2002. This calibration takes the form of 
a matrix with which the Gaussian covariance is multiplied by, to 
obtain the non-Gaussian covariance. This letter is organised as fol- 
lows. The Section 2 provides the notation relevant for this work, 
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and the theoretical description of the Gaussian covariance. Section 
3 describes the ray-tracing simulations and Section 4 shows our 
results. In Section 5 we show their impact on current and future 
contiguous weak lensing surveys. We conclude by discussing the 
limitation of our approach and the work that remains to be done in 
order to achieve percent level accuracy in the non-linear covariance 
estimate. 



2 COSMIC SHEAR AND COVARIANCE 

We follow the notation of Schneider et al. 1998. The power spec- 
trum P K (k) of the convergence ft is given by 
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where /k(w) is the comoving angular diameter distance out to 
a distance w (wh is the horizon distance), and n(w(z)) is the 
redshift distribution of the sources. Ps,D{k) is the 3-dimension 
non-linear mass power spectrum (Peacock & Dodds ~1996l 
ISmith et al. 2003b . and k is the 2-dimension wave vector perpen- 
dicular to the line-of-sight. We are interested in the non-Gaussian 
covariance of the two-point shear correlation function, be- 
cause it can be easily transposed to other two-point statistics 
( Schneider et al. 2002 1 by a suitable integration in fc-space. The 
shear correlation function measured at angular scale i? can be split 
into two components, £±, where 
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and Jo, 4 is a Bessel function of the first kind, of zeroth order for 
and of fourth order for The covariance matrix Cov(£+; $2) 
of the total shear correlation function can be written as a sum of 
three different parts: 

Coir(e+;0i,0a) = <£+(0i)£+OM> = 
D8k{"&i — #2) + q++ + < 4" 1 order correlations > (3) 

The first term is the diagonal statistical noise, depending on the 
intrinsic ellipticity variance, <j e , the total area of the survey, A, and 
the density of galaxies, n. In practical units gives: 
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where A8 is the bin size used for the sampling of the correla- 
tion function. The second term represents the coupling between the 
noise and two point shear correlation function: 

(5) 
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and it can easily be calculated using a prediction for non-linear 
shear power spectrum (Peacock & Dodds 1996 , Smith et al. 2003 ). 
The third term requires the knowledge of the fourth order shear 
correlation function as a function of scale. If we assume Gaus- 
sian statistics, it can be expressed as a sum of two terms 
dSchneider et al. 2002) : 
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dip2(,-(\lp b \) [cos4<£ a cos 4ip b + sin4<fi a sin 4ip b ] 



and tp a , (pt are the polar angles of lp a , lj) b , respectively, cos 4<p a = 
1 - 8^1 ^2 1 1 -0 J 4 , sin A Va = ^ai ^2 ( Vvi - A ) 1 1 V> J 4 , and 
the analogous expressions for ipt,. 

In this paper we are interested in the last term of eq. l[3}- At 
large scales we know that we can use the Gaussian approxima- 
tion and write it as the sum of r + o and r+i. At small scales the 
Gaussian statistics break down and this term cannot be calculated 
with semi-analytical techniques. The rest of the paper discusses our 
technique to calibrate the Gaussian prediction of this quantity in 
order to fit the non-Gaussian value measured in ray-tracing simu- 
lations. Therefore using ray-tracing simulations, we will measure 
the covariance of £+, Cov measU red (£+ ; $1, $2), assuming a e = 0, 
so q++ = and D — and we will define #2), the ratio 

between the measured covariance matrix and Gaussian expectation 
for the covariance matrix: 
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3 DESCRIPTION OF THE SIMULATIONS 

We performed 16 particle in mesh (PM) dark matter simulations 
to cover a light cone of angular size 7x7 degrees, from redshift 
z = to z ~ 3, using the tiling technique proposed by White & 
Hu (2000) and explained in Appendix B of Hamana et al. (2002). 
We used 7 simulations of size 200 Mpc, 4 of size 400 Mpc, 3 of 
size 600 Mpc and 2 of size 800 Mpc. Each Af-body experiment 
involved 256 3 particles in a grid of size 1024 3 to compute the 
forces. The cosmology is a standard ACDM model with = 0.3, 
Obaryons = 0.04, A = 0.7 and H = 70 km/s/Mpc, closed to 
the concordance model (Spergel et al. 2006), with a slightly higher 
value for the normalisation as = 1. Combining the simulation data 
in different ways, we generated 64 different, albeit not fully inde- 
pendent (see below), light cones. Each of them is divided in 64 suc- 
cessive redshift planes separated from each other by 100 Mpc. The 
ray-tracing method is described in Hamana et al. (2002). The spa- 
tial resolution of our simulations translates in an angular resolution 
of the order of 6 ~ 0.5 arcmin for z > 0.2. Given the limitations 
of the PM technique, discreteness effects can be significant at red- 
shift z > 1.5 (due to transcients). Nevertheless, our measurements 
are reliable at scales larger than the mean interparticle distance, i.e. 
8 ~ 2 arcmin. and we expect they can still used with high confi- 
dence level down to 6 ~ 1 arcmin. 

The size, S, of our light cones matches closely that of the sim- 
ulations, so using the dispersion among them to compute the co- 
variance matrix would certainly underestimate its amplitude, even 
at small angular scales. Fluctuations at scales larger than the simu- 
lation box size are also missed with these realisations. Furthermore, 
they are not strictly independent, since they just combine in differ- 
ent ways the 20 simulations. For these two reasons, in the case of 
A — S the value of ^(i^i, $2) on small scales would be always 
underestimated, as compared to the cases A 7^ S, and would not 
converge to unity at large scales. In order to minimise these limita- 
tions and still have a fair estimate of the covariance matrix on the 
estimator used here, it is thus wise to always keep the angular size 
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of the survey A to a small fraction of S. In practice, we divide S 
in 4, 9 and 16 adjacent subsamples, leading to assumed values of 
A ~ 12, 5.4 and 3.1 square degrees and 256, 576, 1024 realisa- 
tions respectively, in total. Note that the choice of A is made such 
that the largest angular scale considered, 6 — 20 arcmin, remains 
small compared to \f~A. We finally choose A — 5.44 deg 2 . 



4 DESCRIPTION OF THE MATRIX CALIBRATION 

We measure .F(i?i,i?2) according to eq. as follows. The term 
Cav me asu*ed(£+;#i! #2) is given by ((£+ - (£+)) 2 }, where £+ is 
measured in each realisation of the survey of size A — 5.44 deg 2 , 
while the average (• ■ ■} is performed over all the realisations. The 
term CovGaussian(£+; $2) is calculated by measuring and 
£_ in the 64 largest samples S of area A — 49 deg 2 , and inte- 
grating numerically eqs. ((6). This ensures that the numerator and 
denominator in eq. are self-consistently defined. It is worth 
noticing that for all cases with A < S the asymptotic behavior 
of .F(#i,i?2) does not converge to unity. It indeed seems to be 
even worse than for the case A = S. This is a well-known effect 
which occurs when the scales become comparable to the size of the 
survey (Peebles 1974). The result is that for those scales the mea- 
sured shear correlation is biased to lower values. Therefore, at small 
scales, the measured cosmic variance Cov meaB ured (£+ ; $1, $2), 
when A < S , is more biased low and decreases faster when 
the scale increases than for the case CovGaussian(£+; #1, #2) and 
A = S. The final result is that the ratio J-($i , #2) becomes smaller 
than unity. Note that in practice, for numerical reasons we have 
to use A = S to compute CovGaussian(£+; , #i 1 i?2) using[6] We 
do not expect this has any impact on our results, within the level 
of accuracy we can achieve from this set of simulations, provided 
we rescale the covariance matrix only in the inner part. The left 
panel of FigQ]shows the diagonal elements , #2) for different 
source redshifts. For a source redshift z s ~ 1 the calibration factor 
is ~ 10 at $1 = 1?2 = 1 arcmin, implying that the cosmic vari- 
ance has been widely underestimated in previous lensing surveys 
at scales below ~ 10 arcminutes. The correction factor is larger for 
lower source redshifts. The transition scale # c , which defines the 
angular scale transition from Gaussian and non-Gaussian covari- 
ance, is redshift dependent because the non-linear regime starts at 
larger scales for nearby structures. Therefore, the calibration ma- 
trix must be parameterized with an explicit redshift dependence. 
We choose a generic power law behavior, as suggested by the left 
panel of Fig.fj] to parameterize .F(l?i, #2): 
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The two panels in Fig. [2] show a and j3 as measured in the 
ray-tracing simulations at nine different source redshifts z s — 
[0.4, 0.5, 0.8, 1.0, 1.2, 1.4, 1.5, 2.0, 3.0]. These measurements are 
well fit by the following redshift dependent functions: 
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For q, we find (01,03,03) = (16.90,0.95,-2.19), and for 0, 
(61,62,63,64) = (1.62,-0.68,-0.68,-0.03) in the samples 
with angular size A — 5.4 deg 2 . The fit is performed on scales 
below 10 arcminutes, which allows us to define also the transition 
angle # c as the scale where the fitted function crosses the Gaussian 
covariance. The third panel of Fig.f2]represents the measurement of 



i? c . Using the same functional form as for a, namely $ c = -i- +t 3 , 
We find the best fit values (h, t 3) ts) = (8.07, 0.95, 1.65). Z 

Since the normalisation of our simulations is high (<T8 = 1), 
we expect $ c to be slightly overestimated. Several other sources of 
uncertainty in our measurements may also spoil the estimate of the 
covariance. In particular, as previously anticipated, there is also a 
"cosmic error" and a "cosmic bias" that affect our measurements 
(e.g. Szapudi & Colombi 1996), which are difficult to estimate. 
Fortunately, such a cosmic bias/error is expected to increase with 
the survey size A. According to eq. (|6), the covariance scales as 
oc 1/A, so T should in fact be independent of A, which allows one 
to use our parametrisation of T for any (reasonable) angular sur- 
vey size. This property can also be used to check the convergence 
between our realisations of various survey sizes as illustrated by 
right panel of Fig. Q] Surveys with areas A = 3.1, 5.4 and 12 
square degrees agree with each other, but there is a problem with 
A = 49 deg 2 , where T is biased low. In the latter case, this is not 
surprising since the light cone size is comparable to the simulations 
size, as discussed in § [3] The convergence between other values of 
A suggest that the cosmic bias/error on T measured in these sam- 
ples is small, i.e. the full set S from which they are extracted, is a 
fair enough sample. We check this by dividing our A = 12 deg 2 set 
of 256 realisations into 4 subsamples of 64 realisations, and mea- 
sured T in each of the subsamples. The dispersion between these 4 
subsamples is of the order of 10% - 20%, which gives a rough idea 
of the accuracy of our estimate of F (■&,&), in agreement as well 
with the convergence between the measurements observed on right 
panel of Fig.[j]for A sC 12 deg 2 . 

While our choice of parametrisation eq. l[8j is globally accu- 
rate to ~ 20% along the diagonal of the matrix J-(fti , #2), it be- 
comes less accurate for very different #1 and #2. One should note 
that the lack of accuracy in the off diagonal components is not criti- 
cal because the cross-correlation coefficient is ,$0.1 in this region. 



5 IMPACT OF NON-GAUSSIANITY ON CURRENT AND 
FUTURE SURVEYS 

Finally, we compare the amplitude of statistical and cosmic vari- 
ance at small scales for a range of contiguous surveys such as 
GEMS ( |Heymans et al. 2005|>, COSMO S (Massey et al. in prep. 
), CFHTLS Wide jHoekstra et al. 20061 and two different versions 
of SNAP (Refregier et al. 2004) whose characteristics are shown in 
table Q] The statistical noise is computed using eq. © , assum- 
ing a bin size A9 = 0.1. Note that the statistical noise differs 
if the bin size used to measure the correlation function is differ- 
ent. In addition we choose a e — 0.4 for ground-based surveys and 
(j e = 0.3 for space-based surveys. Fig. fS] shows that by dropping 
the Gaussian approximation the total noise changes at small scales. 
The changing due to the non-Gaussian correction depends on the 
relative amplitude of the three different contributions to the total 
variance, namely, the shot noise, the sampling noise and the cou- 
pling term. For "low density" surveys, such as the CFHTLS Wide, 
the impact of the non-Gaussian correction is smaller as compared to 
the one expected for the low noise space based surveys, where the 
cosmic variance far exceeds the statistical noise. It is worth notic- 
ing our results are obtained for a higher trg value than Spergel et 
al. 2006 (as — 0.75) and are likely to be slightly different for this 
model. A more extensive analysis of simulations made with differ- 
ent cosmologies would be necessary to accurately predict the am- 
plitudes of the non-Gaussianity corrections to the cosmic variance. 
Unfortunately, only a small set of AC DM ray tracing simulations 
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Figure 1. Left panel: Diagonal elements of the matrix , $2) for different source redshift planes. Right panel: Diagonal elements !F{fti , $2) for different 
survey sizes and z s = 1. The black solid line represents the best-fit of , #2 ) using eq. {8). Error bars are computed using bootstrap with 1000 realisations. 
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Figure 2. The plots shows the measured parameters a, (3 and $ c of the calibration matrix (see eq. [8) as a function of the redshift. Error bars are computed 
using bootstrap with 1000 realisations. The solid line shows the best fit from eq. j9}- 



Table 1. Main Characteristics of surveys used in Fig. [5] 



Name 


A (deg 2 ) 


n 


< Z S > 


GEMS 


0.25 


65 


1 


COSMOS 


1.6 


80 


1.2 


CFHTLS Wide 


50 


15 


0.8 


SNAP deep 


15 


300 


1.1 


SNAP wide 


260 


120 


1.2 



with as = 0.8 is available. This set of simulations, whose char- 
acteristics are given in Heymans et al. 2006, is composed of two 
redshift planes each containing 12 simulations of 25 deg 2 which 
is not enough to find a recalibration fitting formula. Nevertheless, 
Fig|4]shows that even for a A.CDM model with erg = 0.8 the cos- 
mic variance has been widely underestimated. FigOalso shows that 
using a rescaling obtained from rj8 = 1.0 gives results which are 
in good agreement with the ones obtained for as = 0.8 for low 



redshift surveys and slightly overestimates the cosmic variance as 
the depth increases. These simulations were also used to confirm 
the validity of our statements regarding the behavior of the ratio 
J r (i?2, ^2) and the change of the size of A used for the recalibra- 
tion. 



6 DISCUSSION AND CONCLUSION 

We have shown that the non-Gaussian contribution to the covari- 
ance in two-point shear statistics cannot be neglected at small an- 
gular scales. Using ray-tracing simulations we have calibrated the 
non-Gaussian covariance with respect to the Gaussian covariance 
as calculated in Schneider et al. 2002. We have derived a calibra- 
tion matrix which can be used as a first approximation for cosmo- 
logical parameter measurements in current lensing surveys and for 
parameter forecasting. 

We found that the correction coefficient could be as high 10 
at 1 arcminute for a source redshift of 1, and 30 for source redshift 
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Figure 3. The total, statistical plus cosmic variance noise for each survey 
of table [T] The noise including the non-Gaussian correction (thick lines) is 
compared on scales i? 10 arcmin, with the noise expected in the case of 
Gaussian statistics (thin lines). Dropping the Gaussian assumption increases 
noise on small scales. The impact of the non-Gaussian correction for the 
CFHTLS Wide is small; the statistical noise oc 1 /n 2 and the coupling term 
q++ oc 1/n still dominates at small scales for such a density. These same 
terms become negligible for the space based surveys whose density is much 
higher. 
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Figure 4. Diagonal elements of the matrix J r { r &i,'&2) for planes with 
z s ~ 0.5 (blue triangles) and z s ~ 1 (black diamonds) obtained for the set 
of simulations ACDM with ag = 0.8. Error bars are obtained using boot- 
strap with 1000 realisations. For each of the two measurements we compare 
the obtained value with the predicted value calculated from AC DM simu- 
lations with erg = 1., marked with black and blue dotted lines. 



of z s — 0.5. The transition between Gaussian and non-Gaussian 
covariance occurs around 10 arcminutes for z s = 1 and 20 ar- 
cminutes for z s = 0.5. Our work shows that it is important to 
include this non-Gaussian contribution to the shear estimated er- 
rors, and that sub-arcminute resolution ray-tracing simulations are 
very useful for this purpose. Although this source of error has been 
neglected in previous lensing analysis, we note that it should not 
strongly impact the measurement of ag for surveys using the shear 
signal measured above the transition scale # c , where the Gaussian 
covariance is a reasonable assumption. However, it will signifi- 
cantly affect the joined Q]\i-ag constraints, since the degeneracy 
breaking between these two parameters is based on a the relative 
amplitude of the shear correlation signal between small and large 
scales (Jain & Seljak 19971. An increased error at small scale, as 
shown here, will make the degeneracy more difficult to break. 

Extension of this work via a thorough analysis of the non- 
Gaussian covariance based on numerical simulations include shear 
error calibration with broad redshift distribution (tomography), dif- 
ferent two-points statistics and the dependence of the non-Gaussian 
correction with a varying cosmology. In particular we expect a non- 
trivial dependence of the calibration matrix with ag, since, for a 
fixed angular scale, non-linear structures form earlier for higher ag . 
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